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1.  Introduction 

Assume  two  parallel  line  segments  of  indefinite  length,  one  of  which  is  fixed 
and  the  other  of  which  is  movable,  moving  in  a  linear  direction,  parallel  to  the 
fixed  line.  Along  the  fixed  line  there  are  special  points  equally  spaced  which  we 
shall  denote  as  positions.  Similarly,  on  the  moving  line,  there  exist  equally  spaced 
points  which  we  denote  as  sites.  A  site  may  be  in  either  of  two  states  which  we 
refer  to  as  vacant  or  filled.  After  a  site’s  initial  state,  which  is  determined,  it  can 
change  states  only  after  interactions  with  the  positions  along  the  fixed  line  in  the 
following  manner.  A  site  which  is  vacant  can  become  filled  at  certain  positions. 
Call  these  “load”  positions.  Similarly,  a  site  which  is  filled  can  become  vacant 
at  the  remaining  positions.  We  denote  these  as  “release”  positions.  We  assume 
an  arbitrary  fixed  starting  position  on  the  fixed  line  with  the  positions  numbered 
consecutively  beyond  the  starting  position  and,  further,  the  release  and  load 
positions  alternate  so  that  the  odd  numbered  positions  are  release  positions  and 
the  even  numbered  ones  are  load  positions.  Then,  under  certain  assumptions 
stated  in  the  next  section,  the  question  posed  is  the  probability  a  site  will  be  filled 
(or  vacant)  after  a  transit  of  n  positions.  The  model  will  then  be  extended  to  the 
case  where  a  site  can  be  in  any  one  of  (m  +  1)  states  and  the  analogous  question 
posed  is  the  probability  the  site  will  be  empty,  filled  or  in  any  arbitrary  state  j 
after  a  transit  of  n  positions. 

The  model  given  above  is  related  to  the  following  theory  concerning  the  mech¬ 
anism  of  muscular  contraction  as  discussed  by  Podolsky  [1].  A  muscle  fibril,  as 
seen  under  the  electron  microscope,  consists  of  alternating  thick  and  thin  fila¬ 
ments.  It  is  assumed  that  sites  exist  along  the  two  kinds  of  filaments  at  which 
certain  chemical  interactions  occur  at  the  molecular  level.  The  sites  on  the  thin 
filaments  are  capable  of  binding  certain  molecules  and  when  a  site  containing 
the  molecule  approaches  a  site  on  the  thick  filament  an  interaction  (splitting 
off  of  the  molecule)  may  or  may  not  take  place. 

Observations  on  the  living  muscle  suggest  that  during  shortening  there  is 
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relative  motion  of  the  two  filaments.  In  our  model,  we  have  introduced  an  ideal¬ 
ization  by  assuming  that  only  the  thin  filament  moves  (the  moving  line)  and 
that  the  thick  filament  is  rigid  (the  fixed  line).  Further,  Podolsky  assumes  that 
the  positions  on  the  thick  filament  are  only  release  positions  and  that  the  site 
on  the  thin  filament  picks  up  (becomes  filled)  a  molecule  from  the  medium 
between  the  two  filaments  at  any  point  between  two  positions.  For  simplicity, 
we  have  placed  this  binding  action  at  a  fixed  point  between  two  positions  and 
also  call  this  point  a  position. 

2.  Assumptions  and  notation 

Letting  “event”  refer  to  either  the  occurrence  or  nonoccurrence  of  an  inter¬ 
action,  we  make  the  following  assumptions  with  regard  to  the  process: 

(a)  an  event  between  a  site  and  a  position  is  independent  of  any  previous 
events  involving  that  site  and  position, 

(b)  only  one  event  can  occur  between  a  specific  site  and  a  specific  position, 

(c)  the  moving  line  moves  with  uniform  velocity  past  the  fixed  line, 

(d)  a  filled  site  becomes  vacant  at  a  release  position  with  a  constant  probability 
which  is  the  same  for  all  release  positions ;  similarly,  a  vacant  site  becomes  filled 
at  a  load  position  with  a  constant  probability  which  is  the  same  for  all 
load  positions. 

Assumptions  (c)  and  (d)  deviate  from  the  real  situation  in  the  muscular  con¬ 
traction  problem  in  that  during  shortening  a  site  on  the  thin  filament  moves  with 
a  varying  velocity.  The  probability  of  an  interaction  between  a  site  and  a  position 
will  clearly  be  greatly  affected  by  the  speed  with  which  a  site  moves  past  a 
position:  the  slower  the  speed  the  greater  the  probability  of  an  interaction.  Hence, 
a  solution  based  on  assumption  (d)  can  only  be  a  first  approximation  to  the 
question  posed  in  the  muscle  problem.  A  more  realistic  assumption  would  postu¬ 
late  the  probability  of  an  interaction  between  site  and  position  as  a  function  of 
speed  or  time. 

We  let  S  represent  the  state  of  a  site,  *S0  being  the  state  of  the  site  when  it  is 
vacant  and  Si  when  it  is  filled.  If 

(1)  (S,^S,)|n,  =  0, 1;  »  =  1,  2,  •••,(2r  — l),2r,  ••• 

represents  the  transition  of  a  site  from  state  i  to  j  at  the  nth  position  and  recalling 
that  an  odd  position  is  a  release  position  and  an  even  position  is  a  load  position 
we  let 


a  =  -P{(Si  -S„)|2r-  1}, 
0-PftS,  — S0|2r-1}, 
(2)  0  =  P{(S„-^S,)|2r}, 

0  =  P{(S,— »&,)|2r}, 

P„  =  P{(S  =  SOM, 


1  —  a  =  P{(Si  — *  S0|2r  —  1}, 
1  -  P{(S„  -»  S,)|2r  -  1}, 
1  —  0  =  P{(So  — >  So)|2r}, 

1  =  P{(S,  — >  S0|2r}, 

1  -  P.  =  P{(S  =  S„)|n}. 
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We  seek  a  solution  for  P„,  the  probability  the  site  will  be  filled  at  the  end  of 
the  nth  position. 


3.  Solution  through  difference  equations 

The  branching  process  can  be  illustrated  as  shown  in  figure  1 . 

(2r  —  l)st  Position  2rth  Position 

(Release) 


Figure  1 

We  write  at  once  the  two  joint  difference  equations 

(3)  P2r  —  P  2r— 1  +  /3(1  —  P  2r— l)  =  0  +  (1  —  0)Pir-l 

(4)  P  2r— 1  =  (1  Oi)P  2r— 2 

from  which  we  obtain  the  following  second  order  difference  equations 

(5)  P2r  =  0  +  (1  -  o)(l  -  0) P 2r— 2 

(6)  Pi  r-i  =  /3(1  —  a)  +  (1  —  a)(l  —  /3)P2r_3. 

The  standard  solution  of  (5)  and  (6)  gives  respectively 

(7)  P„  =  [d  -  «)(1  -  0)]'P.  +  ft 

and 

(8)  P*_,  =  [(1  -  «)(1  -  | B)]«i>,  +  «1  -  oc)  1  1~I<1(1~J,)a()1(1~_y 

=  (1  -  «)[(1  -  a)(l  +  -a)  -  ~_y  ’ 

r=  1,2,  •••  . 

Thus,  if  the  site  starts  in  state  Si,  the  boundary  condition  is  P0  =  1  and 

(9)  P.  =  [(1  -  a)(l  -  «]"'*  +  0  - 

n  =  0,  2,  4,  •  •  • 


r  =  1,2, 


or 
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(10)  P„  =  (1 -<*)[(!- a)(l 


+  jg  -  a)  '—.ULr  «)fl  -  m 

f  PU  OL)  1  _  (1  _  _ 


(1  -  a)(l 

If  the  site  starts  in  state  S0,  so  that  P0  =  0,  we  have 


(»— 1)/2 
0T~’ 


(ii) 


i  -  ra  -a)g  ~^)i"/2 

"  P  1  -  (1  -  a)(l  -  0)  ’ 


7i  —  1, 3,  5,  *  *  * . 

n  =  0,  2,  4, 


or 

(12) 


Pn  =  0(1  —  a) 


i  -  ra  -«)a  -0)i(n-»^2 
i-a-^a-0)  * 


7i  —  1,  3,  5 ,  •  •  • . 


4.  Solution  through  transition  matrices 

It  is  of  interest  to  obtain  the  solution  to  this  problem  involving  only  dichoto¬ 
mous  states  of  a  site  using  transition  matrices.  This  method  can  be  employed 
to  solve  the  more  general  problem. 

From  the  branching  process  given  previously,  we  obtain  the  matrix  of  transi¬ 
tion  probabilities  for  the  state  of  a  site  as  it  passes  an  odd  position  to  be  shown 
in  table  I,  and  the  transition  matrix  as  the  site  passes  an  even  position  shown 
in  table  II. 


TABLE  I 


State  before 

State  after  (2 r 

—  l)st  position 

(2 r  —  l)st  position 

So 

Sr 

So 

1 

0 

Sr 

a 

1  —  a 

TABLE  II 


State  before 

State  after  2rth  position 

2rth  position 

So 

Sr 

So 

1  -  0 

P 

Sr 

0 

1 

Denote  the  transpose  of  the  above  matrices  by  A  and  B  respectively,  that  is, 

<13>  A  =  (o  l-«) 

and 


Let  Vi, n  =  P{(S  =  Sj)\7i},  withj  =  0,  1  and  P'n)  =  (po,n,  Pi.n),  a  row  vector  of 
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probabilities  of  the  two  possible  states  of  a  site  at  the  end  of  the  nth  position. 
Then 


(15) 

P(2r)  = 

B  P  (2r— 1) 

P(2r-1)  = 

A  P  (2r— 2), 

Hence  for  n  even, 

(16) 

P(n)  = 

=  BAP(„_ 

2)  =  (BA)n/2P(0). 

Since 

(17) 

BA  = 

n  -$ 

\  0 

«(1  -  0)  \ 
)) 

we  have 

(18) 

ft:)  -  (’ 

L  -0 

0  1 

o(i  -  p)  y2 

-  a(l  -  P)J 

n  =  0,  2,  4,  •  •  •  , 

where  po.o  and  pi.o  are  the  initial  conditions  of  the  process.  If  the  site  starts 
filled,  po.o  =  0  and  pi.o  =  1;  if  it  starts  empty  po,o  =  1  and  pi,o  =  0. 

For  n  odd,  we  get 


(19) 

which  yields 

(20) 


Poo  =  ABP(n_2)  =  (AB)  (n-1)/2P(1) 

(Po,n\  _  /I  —  0(1  —  a)  a  \(n_1)/2/p0li\ 

\2W  \  /8(l-o)  1  ~  a)  \pi,J 


n  —  1,  3,  5,  •  •  •  , 

where  po,i  =  a  and  pi,i  =  1  —  a  if  the  site  starts  filled,  and  po.i  =  1,  Pi,i  =  0 
if  the  site  starts  empty. 


5.  Generalization 


We  extend  the  problem  to  the  case  where  each  site  can  be  in  any  one  of  m  +  1 
states  /S0,  Si,  •  •  •  ,  Sm.  When  a  site  in  state  Si,  with  i  =  0,  1,  •  •  •  ,  m  on  the 
moving  line  encounters  an  odd  (release)  position  on  the  fixed  line  segment,  any 
number  i  —  j,  where  j  =  0, 1,  ♦•*,«,  of  interactions  can  occur  resulting  in  the 
site  changing  from  state  S » to  S,.  Similarly,  when  a  site  in  state  S t-  with  i  =  0, 
1,  •  •  •  ,  m  encounters  an  even  (load)  position  on  the  fixed  line  segment,  any 
number  j  —  i,  where  j  =  i,  i  +  1,  •  •  •  ,  m,  of  interactions  can  occur  resulting 
in  the  site  changing  from  state  S to  S}.  It  is  assumed  that  the  number  of  inter¬ 
actions  occurring  at  the  juxtaposition  of  a  site  and  a  position  has  a  binomial 
distribution  with  parameter  a  at  a  release  position  and  parameter  p  at  a  load 
position. 

We  now  let 


(21) 


Pi.,  =  P{(S  = 


i  =  0,  1,  •  •  •  ,  m, 
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and 


(22)  P(»)  —  (Po,»>  Pl.n,  *  *  *  >  Pm.nY 

be  an  (ra  +  1)  X  1  column  vector.  In  general,  we  desire  to  determine  any  com¬ 
ponent  of  P  (n>  although  there  may  be  some  particular  interest  associated  with 
Po.n  or  pm,n,  the  probabilities  the  site  is  empty  or  filled,  respectively,  as  the  site 
has  moved  past  n  positions. 

Since  we  assume  the  number  of  interactions  is  a  binomial  variate,  we  may  then 
write  for  the  transition  probability  at  an  odd  position 

(23)  P{(Si  ->  Sj)\2r  -  1}  =  {G  -  ”  a)J  =  J  =  h 

10,  j  >  i, 

i  =  0, 1,  •  •  •  ,  m;  j  =  0, 1,  •  •  •  ,  i 


and  for  the  transition  probability  at  an  even  position, 


(24) 


pm 


Sj)\2r}  = 


(mj  -ty’-'i1  ~  ^)w■,  =  ./-<(£)> 

0, 

i  =  0, 1,  •  •  •  ,  m;j  =  i,  i  +  1, 


j  ^  i 
j  <  i, 
,  w. 


Then  the  transition  matrix  at  an  odd  position  is  as  shown  in  table  III, 


TABLE  III 


So  ■ 

•  Si  . 

•  Si  •• 

•  sm 

So 

1  • 

•  0 

•  0 

•  0 

Si 

&..<(«)  * 

•  bi,  o(a)  • 

•  0 

Sm 

bm.m(oO  * 

■bm.m— j(ot)  • 

■bm.m— »(tt)  * 

•bm.o(a) 

where  0&M-i(«)  =  1-  Denote  this  matrix  by  A',  the  prime  representing  the 
transpose. 

The  transition  matrix  at  an  even  position  is  as  shown  in  table  IV, 
where  =  1*  Denote  this  matrix  by  B'. 

We  then  have 

(25)  P(2r)  =  B  P(2r-1) 
and 

(26)  P(2f-1)  =  A  P  (2r— 2) 
so  that 

(27) 


P(2r)  =  BA  P(2r-2)  =  (BA)2P(2r_4)  =  (BA) ^r)/2P(0). 
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TABLE 

IV 

S0  . 

■■Si  • 

•  •  s ,  . 

••  sm 

-So 

6m,o03)  • 

•  •  bm.iOS)  • 

••  bm.M  • 

■■  bm.m  (0) 

Si 

0  • 

’  '&m-».o(/3)  * 

*  •($) 

sm 

0  • 

,.  o  • 

0 

l 

Similarly, 

(28)  P(2r_1}  =  ABP(2r_8)  =  (AB)2P(2r_6)  =  (AB)<2-»/2P(1)  =  (AB)  ^A  P(0). 
Therefore, 

(29)  P(»)  =  (BA)n/2P(0),  w  =  0,  2,  4,  •  •  •  , 

and 

(30)  P(„)  =  (AB)  <»-”/2A  P(0),  n  =  1,  3,  5,  •  •  •  . 

P(0)  is  the  vector  whose  components  {pli0} ,  with  i  =  0, 1,  •  •  •  ,  m,  represent  the 
initial  conditions.  Of  interest,  are  the  cases  where  pm, o  =  1,  pt-,o  =  0,  i  m,  or 
where  po.o  =  1,  p.,o  =  0,  i  ^  0.  The  former  case  holds  when  the  site  starts  in 
state  Sm,  that  is,  it  is  filled,  the  latter  case  represents  the  situation  where  the 
site  starts  in  state  So,  that  is,  it  is  empty. 

6.  Estimation  of  parameters  in  the  simple  case 

If  a  and  0  are  known  then  clearly  the  probability  that  a  site  will  be  either 
empty  or  filled  at  any  position  n  is  immediately  determined  by  equations 
(9)-(12)  or  (16)  and  (19)  in  the  special  case.  Likewise  the  probability  that  a 
site  is  in  any  state  So  to  Sm  is  determined  by  (29)  and  (30)  in  the  general  case. 
If  a  fixed  number  of  sites  in  known  starting  states  are  observed  at  the  end  of  a 
transit  past  n  positions,  then  the  set  of  proportions,  (p  where  i  =  0, 1,  •  •  •  , 
m,  of  sites  observed  to  be  in  state  i  are  random  variables  having  a  joint  multi¬ 
nomial  distribution  or  a  binomial  distribution  when  fi  =  1.  The  statistical 
problems  might  then  be  in  testing  the  binomial  model. 

Of  greater  interest  is  the  case  where  a  and  $  are  unknown.  We  consider  the 
problem  of  estimating  a  and  /3  only  in  the  special  case  where  the  possible  states 
of  the  site  are  dichotomous  so  that  we  have  solutions  in  closed  form  as  given 
by  equations  (9)  to  (12).  We  illustrate  the  procedure  by  employing  equation  (9), 
similar  techniques  being  applicable  to  the  other  forms. 

In  (9),  let  0  =  (1  —  a)(l  —  0)  to  obtain 


(31) 


1  —  Bnl 2 

Pn  =  0"/2+i8ij-=T, 


n  —  0,  2,  4,  •  •  •  ;  0  <  d  <  1. 


264 


FOURTH  BERKELEY  SYMPOSIUM:  GREENHOUSE 


Now  limn_»oe  Pn  =  /3/(l  —  0).  Let  N  be  large  enough  so  that  we  may  write  to 
a  given  degree  of  accuracy 


(32) 

Then 

(33) 

yielding 

(34) 

and  hence 


Pn  = 


<8 

1  -  0 


Pn  =  0"'2  +  PN(  1  -  0"'2) 


0n/2  _ 


Pn  —  Pn 
1  —  Pn 


Let  u  and  v  be  two  independent  estimates  of  P„  and  Pn  respectively.  Then  an 
estimate  of  0  is 

(36)  t  = 

From  (32),  and  the  definition  of  0 ,  we  get 

(37)  i  =  (1  -  6)v 


and 

(38) 


<$  = 


1 


$ 

1-/5’ 


If  n  =  2,  that  is,  u  is  observed  at  the  second  position  on  the  fixed  line,  these 
estimates  become 


(39) 


e  = 


u  —  V 
1  —  V 


(40) 


A 

/3  = 


1  —  u 

1  —  V 


V 


(41) 


=  (1  -  m)(1  -  t>). 
1  —  2v  +  uv 


The  observables  u  and  v  can  be  obtained  in  the  following  manner.  Consider 
starting  with  m  sites  all  in  state  Si,  since  we  are  using  (9)  which  is  obtained  with 
this  initial  condition.  Then  if  the  states  of  each  site  are  noted  at  a  selected  nth 
position,  u  is  the  proportion  of  sites  observed  to  be  in  state  Si.  Repeat  this  with 
another  set  of  sites  and  now  if  the  state  of  each  of  these  sites  are  noted  at  a  large 
distance  along  the  fixed  line  from  the  starting  position  so  that  N  is  large  enough 
for  the  approximation  (32)  to  take  hold,  v  will  be  the  proportion  of  sites  in  Si 
at  the  Nth.  position.  Thus  u  and  v  are  independent,  binomial  variables  with 
expectations  Pn  and  Pn  and  variances  P„(l  —  P„)/m(n )  and  Pjv(1  —  Pjv)/ra(iV) 
respectively. 
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If  the  number  of  sites  involved  in  the  experiments,  m(ri)  and  m(N),  are  large, 
then  we  may  obtain  an  estimate  of  the  variance  of  6  by  the  usual  Taylor  series 
approximation.  This  gives 


(42) 


V(6)  = 


4 

n20(n-2)(l  ~  v)2 


(1  ~  u)2  -2 
(1  -  v)2 


The  estimates  of  variance  of  &  and  $  are  more  difficult  to  obtain  for  general 
values  of  n  because  of  the  covariance  between  v  and  8  and  between  0  and  0. 
However,  for  n  —  2,  the  estimates  &  and  0  are  functions  only  of  u  and  v,  and 
estimates  of  their  variances  are 


(43) 

and 

(44) 


1 

(1  -  v)2 


(1  -  u)2 

(1  -  V)2 


1 

(1  —  2v  +  uvY 


[(i  -  vy&l  +  (l  -  uyrtl 


7.  Concluding  remarks 

The  previous  discussion  was  based  on  two  major  assumptions:  (a)  the  proba¬ 
bilities  of  the  two  kind  of  interactions,  a  and  0,  were  constant  over  all  positions, 
and  (b)  the  transition  probabilities  were  binomial.  We  could  relax  both  conditions 
and  with  little  change  in  the  development  obtain  a  solution  of  P(„). 

Let 

45)  P{(Sf  ->  S,-)|2r  -  1}  =  Pij[a(2r  -  1)], 

i  ^  j,  i  =  0,  1,  •  •  •  ,  m\  j  =  0,  1,  •  •  •  ,  i,  where  a,  the  probability  of  a  release 
interaction,  is  now  a  function  of  the  position  number  along  the  fixed  line,  and  p 
is  some  discrete  probability  function  whose  parameter  is  a(2r  —  1).  Also  denote 
the  transition  matrix  whose  elements  are  the  transition  probabilities  just  defined 
by  A(2r— i).  Similarly,  we  let 

(46)  P{(S,  -  S,)\2r}  =  qii[fi(2r)l 

i  ^  j,  i  —  0,  1,  •  •  •  ,  m,  j  =  i,i  +  1,  •  •  •  ,  m,  where  /3,  the  probability  of  a  load 
interaction,  also  depends  upon  the  position  number,  and  q  is  a  discrete  probability 
function  with  parameter  0(2 r),  not  necessarily  of  the  same  type  as  p.  Denote 
the  transition  matrix  at  an  even  site  by  B2r.  Then  we  have,  for  example,  for  n 
even,  the  very  general  solution 

(47)  P(„)  =  B(n)A(„_i)B(n_2)A<„_3)  •  •  •  B(2)A(i)P(0). 

Now  with  a  and  0  functions  of  n,  we  may  approximate  the  muscle  problem  more 
closely  since  in  effect  a  and  0  can  be  taken  as  velocity  functions. 
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